Damage spreading for one-dimensional, non-equilibrium models with parity 

conserving phase transitions 
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The damage spreading (DS) transitions of two one-dimensional stochastic cellular automata sug- 
gested by Grassberger (A and B) and the kinetic Ising model of Menyhard (NEKIM) have been 
investigated. These non-equilibrium models exhibit non-directed percolation universality class con- 
tinuous phase transition to absorbing states, exhibit parity conservation (PC) law of kinks and have 
chaotic to non-chaotic DS phase transitions, too. The relation of the critical point and the damage 
spreading point has been explored with numerical simulations. For model B the two transition 
points are well separated and directed percolation universality was found both for spin damage and 
kink damages in spite of the conservation of damage variables modulo 2 in the latter case. For 
model A and NEKIM the two transition points coincide with drastic effects on the damage of spin 
and kink variables showing different time dependent behaviours. While the kink DS transition of 
these two models shows regular PC class universality, the spin damage of them exhibits a discon- 
tinuous phase transition with compact clusters and PC-like spreading exponents. In the latter case 
the static exponents determined by finite size scaling are consistent with that of the spins of the 
NEKIM model at the PC transition point. The generalised hyper-scaling law is satisfied. Detailed 
discussion is given concerning the dependence of DS on initial conditions especially for the A model 
case, where extremely long relaxation time was found. 



I. INTRODUCTION 

While damage spreading (DS) was first introduced in 
biology m it has become an interesting topic in physics as 
well f2|-[|] . The main question is if a damage introduced 
in a dynamical system survives or disappears. To investi- 
gate this the usual technique is to make replica(s) of the 
original system and let them evolving with the same dy- 
namics and external noise. This method has been found 
to be very useful to measure accurately dynamical expo- 
nents of equilibrium systems [El . It has turned out how- 
ever, that the DS properties do depend on the applied dy- 
namics. An example is the case of the two-dimensional 
Ising model with heath-bath algorithm versus Glauber 
dynamics |i] |j. 

To avoid the dependences on dynamics Hinrichsen et 
al. H suggested a definition of "physical" family of DS 
dynamics according to which the active phase may be 
divided to a sub-phase within which DS occurs for ev- 
ery member of the family, another sub-phase where the 
damage heals for every member of the family, and a third 
possible sub-phase where DS is possible for some mem- 
bers and the damage disappears for other members. The 
family of possible DS dynamics is defined such that it is 
to be consistent with the physics of the single replicas 
(symmetries, interaction ranges etc.). 

The universality of continuous DS transitions is an 
other open question. There is a hypothesis raised 
by Grassberger pOl, that damage spreading transitions 



generically belong to the universality class of directed 
percolation (DP) if they are separated from the ordinary 
critical point. This claim is based on the DP hypoth- 
esis applied for the absorbing type DS transition since 
we can consider the difference of the replicas as an other 
dynamical system evolving by a complex rule. Accord- 
ing to the DP hypothesis - conjectured first in the early 
eighties [ [TT||T^ 1 - in the absence of conservation laws ev- 
ery continuous phase transition of a system with scalar 
order parameter and local interactions to a single absorb- 
ing phase would belong the universality class of the DP. 
There are other more complex models like those with 
several absorbing states [|13| and multi-component sys- 
tems |Q, which exhibit DP transition too. During the 
years the DP universality class has been proven to be 
extremely robust. For a long time only a few number 
of exceptions has been found, which don't belong to the 
DP class. These are the parity conservation (PC) models 
and the multiplicative noise systems fl5|| . 

The first examples of the PC models were Grass- 
berger's (A and B) stochastic cellular automata (SCA) 
pf. The kinks ('00' and 'll'-s) of these models exhibit 
mod 2 parity conservation and the absorbing state is dou- 
bly degenerated. Following that a series of models in the 
same universality class have been discovered. In case of 
the branching annihilating random walk with even num- 
ber of offsprings (BAWe) 17 ll| the parity of the "parti- 
cles" is conserved and there is a single absorbing state. In 
the non-equilibrium Ising model with combined spin-flip 
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and spin exchange dynamics (NEKIM) JT^,^] the kinks 
have local parity conserving symmetry as well as the ab- 
sorbing state is symmetrically doubly degenerated. The 
three species monomer- monomer [^lj (3MM) and the in- 
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are multi 



teracting monomer-dimer (IMD) models 
component models with parity conservations and sym- 
metric absorbing phases. 

The common feature of these models that force them 
in a non-DP universality class was first conjectured to be 
the parity conservation (the PC class name comes from 
here). However this had to be refined, because models 
were found with global parity conservation but DP class 
phase transition (see example ref. |22| (ISCA)). Field 
theoretical investigations of the BAW models showed 
that the BAWe parity conservation dynamics in one- 
dimension results in a new non-DP fixed point possessing 
the PC class universality, while for for odd number of off- 
springs (BAWo) the transition is in the DP class p5| , p6]| . 
Furthermore, among the multiple absorbing state models 
one can have DP behaviour and BAWe parity conserva- 
tion dynamics together if the symmetry of the absorbing 
states is broken (see [Q in case of IMD, [^?J in case 
of the NEKIM model and (3lj in case of the General- 
ized Domany-Kinzel SCA (GDK)). This implies that for 
multi-absorbing state models the BAWe parity conser- 
vation is not a sufficient condition to have non-DP uni- 
versality class but the symmetry of the ground state is 
necessary too. 



Summary of PC related models 
The notation "k" refers to kink 
"h" means external field as variable of the model 



"s" refers to spin, 



Model 


Abs. state 


Dynamics 


Univ. 


Rcf. no. 


A-k 


Z2 symm. 


BAWe 


PC 
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B-k 


Zi symm. 


BAWe 


PC 
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PC 
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BAWo-s 
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DP 
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NEKIM-k 


Z2 symm. 


BAWe 


PC 
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Zi symm. 
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PC 
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Zi symm. 


BAWe 


PC 
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,P| 


GDK-s 


Z2 symm. 


BAWe 


PC 
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ISCA-s 


singlet 


Global-PC 


DP 
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IMD-s+h 


Z2 broken 


BAWe 


DP 
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NEKIM-k+h 


Z2 broken 


BAWe 


DP 
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GDK-s+h 


Z2 broken 


BAWe 


DP 
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Grassberger's A and B stochastic cellular automata 
(SCA) pi with synchronous dynamics and the NEKIM 
model [19[. The 'kink' variables of these models possess 
parity conservation and continuous PC class phase tran- 
sition. 

In case of the NEKIM model the '01' and TO' pairs act 
as kinks and follow the basic (BAWe) elementary reac- 
tions : 

• left-right diffusion, 

• X — > 3X reproduction, 

• 2X — » annihilation. 

Since at damage spreading problems we follow the evo- 
lution of two or more replicas, we can consider it as a 
special, multi-component dynamical problem (here with 
multiple absorbing states too). Furthermore, when the 
DS point is inside the active phase there is a passive state 
of the damage variables with fluctuating replicas in the 
background. Therefore the simple DP universality hy- 
pothesis can not be applied here (although this does not 
exclude a DP class phase transition). 

We also studied here the DS properties of spins of these 
systems, which can be regarded as the "dual" variables of 
kinks but they don't obey the parity conservation. The 
combined observation of spin and kink damage variables 
sheds light on the interplay of parity conservation, ab- 
sorbing state symmetry and universality. In this paper 
we shall show by numerical simulations that the univer- 
sality here is determined by not only the dynamics but 
the symmetry of the absorbing state is a necessary con- 
dition again like in cases of ||j,|2^] . 

If the critical point and the DS point coincide inter- 
esting things happen. While the kink damage exponents 
will belong to the PC universality class, in case of the 
spin damage the static exponents determined by finite 
size scaling are in agreement with that of the the pure 
NEKIM model at the PC transition point on the spin 
level . Detailed discussion on this is forthcoming |33) . 

The dependence of the DS results on the initial states 
of the replicas is discussed because for model-A a very 
slow relaxation makes it an important point. 



II. THE DAMAGE SPREADING SIMULATION 
METHODS 



A very recent study has shown that DS transition 
is possible in a one-dimensional non-equilibrium kinetic 
Ising model j3^| too, and the universality class of the 
transition is not always in the DP class. The dynam- 
ics was engineered as the combination of two subrules 
such that it creates Z2 symmetric passive states, the kink 
damage variables follow BAWe parity preserving dynam- 
ics and a PC universality class DS transition emerges. 

In this work we have investigated the damage spread- 
ing behaviour of some one-dimensional PC models: 



The time dependent simulation is a well established 
method to locate critical points and to measure dynam- 
ical critical exponents at the same time [ p8| . Here we 
applied it for kink and spin damage variables for system 
sizes L = 4096 — 16384 with periodic boundary condi- 
tions. A single spin-flip difference is introduced between 
two identical replicas at the beginning of each simulation 
runs. The difference of spin and kink variables is mea- 
sured during a time evolution with identical rules and 
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random numbers for both replica. The maximum num- 
ber of simulation steps was chosen to be tuAX = L/2, 
and so the damage variables can not reach the boundaries 
and one can avoid finite size effects of them. However 
simulating near the critical point causes long transients 
hindering to see the true scaling behaviour within reach- 
able times. 

The role of initial states of the replicas is not discussed 
in DS simulation literature. If the DS transition point is 
not in the neighbourhood of a critical point an exponen- 
tially quick transient to the steady state is expected, but 
if they coincide - as in case of the Grassberger A model 
- the evolution to steady state slows down to power law 
time dependence and we can expect finite time effects. 
First random states have been chosen with equal and 
uniform distribution of 0-s and 1-s. In case of model A 
SCA this resulted in very confusing results. Then we 
investigated the effects from starting with steady state 
configurations i.e. replicas were driven to steady state 
before the DS measurements. 

The quantities characterising damage evolution show 
power like behaviour in the t — > oo limit at the damage 
spreading point (j>d) separating chaotic and non-chaotic 
phases. The Hamming distance will be the order param- 
eter of this paper: 
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where s(i) may denote now spin or kink variables. Kink 
variables for these models are the '00' and '11' pairs in 
case of the A,B SCA and '01' and TO' pairs for the 
NEKIM model. If there is a phase transition point, the 
Hamming distance behaves in a power law manner at 
that point: 
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(2) 



Similarly the survival probability of damage variables be- 
haves as: 
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and the average mean square distance of damage spread- 
ing from the center scales as: 
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The evolution runs were averaged over N s independent 
runs for each different value of p in the vicinity of pd ( 
but for R 2 (t) only over the surviving runs ). 

To estimate the critical exponents and the transition 
points together we determined the local slopes of the scal- 
ing variables. For example in case of the survival proba- 
bility: 



S p (t) 



In [P a (t)/P 3 (t/m)] 
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(5) 



and we have used m = 4. In the case of power-law 
behaviour we should see a horizontal straight line as 
1/t — ► 0, whenp = pd- The off-critical curves should pos- 
sess curvature. Curves corresponding to p > pd should 
veer upward, curves with p < pd should veer downward. 

The damage spreading measurement of the order pa- 
rameter time scaling can be very effectively parallelised 
in a multi-spin code manner |Io[| , since one needs only one 
random number for each site of the different replicas and 
so one can follow the evolution of N r = 32 * 31 replicas 
in a simple 32-bit computer vector of length L. However 
this method is not applicable to measure the survival 
probability scaling and z, since the healing of differences 
among all the N r replicas takes a very long time and one 
can not introduce a single initial damage for each pair at 
the center of the lattice. For the simulation of survival 
probability a very effective code has been implemented 
for a special, associative string processor ]34| . 

To determine static exponents finite-size scaling (FSS) 
simulations were performed as well. As it was shown 
by Aukrust et. al. §| and @, FSS is applicable to 
continuous, non-equilibrium phase transitions. At the 
critical point the order parameter steady state density 
(D) and the fluctuation x = L d (< D 2 > - < D > 2 ) 
scale with the system size as : 
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where v± is the correlation length exponent in the space 
direction : 



£(p) <x \P ~Pc 



(8) 



(3 is the order parameter exponent in the steady state: 
D{p)^\p-p c f , (9) 



and 7 describes the fluctuation of it: 
X(p) x \P~Pc\~ 1 - 



(10) 



Simulations were done in one dimension for lattice sizes : 
L = 64, 128, 2048. The necessary time steps to reach 
steady state were determined experimentally. The time 
evolution of the concentration was plotted, and the nec- 
essary time steps was fixed for a given L such as 200 — 500 
time values following the level off. Averaging was done 
for these 200 — 500 values times the number of surviv- 
ing samples (500). The pd values were taken from the 
time-dependent MC calculations. 

The dynamic exponent Z = v\\/v± can be determined 
from the FSS of the characteristic time r(p,L). In this 
study we measured the time necessary to reach half of the 
steady state concentration starting from a single damage 
state. The characteristic time obeys the finite size scaling 
law : 
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-(p,L)<xL z h((p-p c )LV^) ; 



(11) 



where Z = V\\jv^. For this measurement we used the 
same damage concentration time evolutions as in case of 
the static runs above. 



III. THE GRASSBERGER B MODEL 

A (BAWe) parity conserving dynamics can be realized 
on the kinks (or 'particles') of the following SCA (we 
show the configurations at t — 1 and the probability of 
getting '1' at time t): 

t-1: 100 001 101 110 011 111 000 010 
t: lllppOOO 

The '00' and '11' pairs are the simplest kinks of the 
model. The time evolution pattern for p = is a reg- 
ular chess-board in 1 + 1 dimension (Rule-50 with double 
degeneration), i.e. , the absorbing states are period-two 
antiferromagnetic. For p < p c (— 0.539(1)) kinks disap- 
pear exponentially, while for p > p c they survive with a 
finite concentration. In the p = 1 limit we get the de- 
terministic Rule-122, which is known to be chaotic. So 
there is damage spreading phase transition besides the 
absorbing phase transition of PC universality. 



A. Kink damage results 

First the simulations were started from two replicas of 
lattices with identical random initial states but with a 
2-kink initial difference. The parity of the lattice forces 
even or (odd) number of initial kinks, therefore it is not 
possible to create odd numbered kink differences. The 
parity of kinks is conserved. The parity of kink differ- 
ences (even) is conserved too. 
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FIG. 1. Local slopes of the Hamming distance (rj) in the B 
model, for p = 0.62, 0.63, 0.632, 0.634, 0.65 (curves from bot- 
tom to top). Statistical averaging was done over 10000 sam- 
ples. 



Still we see a DP like universality of the damage vari- 
ables (Fig. 0,§). The location of the damage spreading 
point (pd — 0.632(1)) is far from the PC critical point 
(p c = 0.539(1)), therefore the active phase is divided to 
a chaotic and non-chaotic sub-phase similarly to the case 
of the Domany-Kinzel SCA [||33|LCj]. 
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FIG. 2. Local slopes of the damage survival probability (S) 
in the B model, for p = 0.62, 0.63, 0.632, 0.634, 0.65 (curves 
from bottom to top). Statistical averaging was done over 
10000 samples. 

The replicas at the DS point are in fluctuating states 
therefore they have not the " chess-board" double degen- 
eration as in case of the PC critical point. The emergence 
of the DP exponents in spite of the mod 2 conservation 
of kink damage variables is similar to what was found 
numerically and analytically for PC models, when the 
external H field destroyed the symmetry of the absorb- 
ing state p4|j27|j3l|]. This suggests that the BAWe parity 
conservation rule is not a sufficient condition for having 
non-DP universality. 

The same results has been found if we started the repli- 
cas from steady states with a single spin-flip initial dif- 
ference. 



B. Spin damage 

Following the damage (difference) of the spins, instead 
of the kinks, we obtained the same DP like results. The 
universality was insensitive to the parity of the damage 
variables, it is DP for both cases. 



C. Finite size scaling for both cases 

Finite size scaling simulations were performed at the 
DS transition point (pd ~ 0.633) for system sizes L — 
64, 128, ...1024. The necessary time steps to reach steady 
state were t = 40000,80000, ... respectively. The results 
can be seen on Figure |^. As on can observe both the 
kink and spin damage concentration show a scaling with 



4 



—f3/ l/ ± — —0.25, while the fluctuations have a slope 
j/i/± = 0.5, all agreeing with the corresponding expo- 
nents in the DP universality class. For the critical dy- 
namical exponent Z = v\\/v±_ DP-like scaling has been 
found for both the spin and kink damage cases. Fitting 
can be done with i 1 - 5798 (DP value) in both cases. 
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FIG. 3. Finite size scaling results for the spin and kink 
damage for the B model. The diamonds correspond to the 
kink damage concentration, the squares to the fluctuations of 
it. The crosses correspond to the spin damage concentration, 
the "x"-s to the fluctuations of it. Triangles and stars denote 
characteristic times r of the spin and the kink damage cases. 
Averaging was done over 500 surviving samples. 



IV. THE GRASSBERGER A MODEL 

Another very similar model exhibiting parity conser- 
vation of kinks is the Grassberger A stochastic cellular 
automaton: 

t-1: 100 001 101 110 Oil 111 000 010 
t: 1 1 1-p 1-p 1 

The time evolution pattern in 1 + 1 dimension, for small 
p evolves towards a stripe-like ordered steady state (with 
double degeneration), while for p > p c (0. 1245(5)) the 
kinks (the '00' and '11' pairs survive. For p — we 
have the Rule-94, class 1 CA, while the p = 1 limit is 
the chaotic Rule-22 deterministic CA. Therefore we can 
expect a damage spreading phase transition between p = 
p c and p — 1 (damages always heal or survive in the 
ordered steady states). 



A. Kink damage results 

First two replicas of lattices of the same random initial 
distributions but a single spin-flips difference condition 
were followed. As the local slope figure (Fig. Eh of the 
D(t) shows the DS transition point (pd — 0.133(1)) seems 
to be slightly off the critical point (p c — 0.1242(5)). 
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FIG. 4. Local slopes of the Hamming distance {rf) in the 
A model, for p — 0.130,0.132,0.134 (curves from bottom to 
top). The simulations were started from random initial state. 
Statistical averaging was done over 6 x 10 5 samples. 

One can read off r\ = 0.31(1), which is close to the 
DP universality class value (j\bp = 0.314(3)), but for the 
survival probability we got nearly zero exponent (Fig. |^). 
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FIG. 5. Local slopes of the damage survival probability 
(5) in the A model, for p = 0.130,0.132,0.134 (curves from 
bottom to top). The simulations were started from random 
initial state. Statistical averaging was done over 10 4 samples. 

The survival probability scaling with i5 ~ contradicts 
the DP scaling and one may speculate, that we can see a 
finite time effect. We extended the same time dependent 
simulations up to t max = 28000 for certain p values, but 
there were no sign of change in the above results. 

To check the transients the evolution of the kink con- 
centration starting from a disordered state have been fol- 
lowed on a L = 8192 lattice. As the Figure ^ shows there 
is a very long relaxation in this model, and the steady 
state has been reached following 2 x 10 6 MC time steps 
only. 

Therefore time dependent simulations from steady 
state initial conditions have been performed. The ini- 
tial states now were chosen to be the outcomes of runs 
following 5 x 10 6 time steps for different p-s, with the 
usual single spin-flip difference. Now we can see dramat- 
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ical changes. First, the DS point moves to the critical 
point (p d = 0.1242(l))(see Fig. @). 
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FIG. 6. The evolution of the kink concentration in the 
A model (L = 8192) started from random initial state at 
p = 0.124. 



-0.26 



-0.28 




-0.32 

0.002 0.004 0.006 0.008 0.01 
1/t 

FIG. 8. Local slopes of the survival probability (8) in the A 
model, for p = 0.123, 0.124, 0.1245, 0.125 (curves from bottom 
to top) in case of steady state initial condition. Averaging was 
done over 50000 independent samples. 



B. Spin damage results 
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FIG. 7. Local slopes of the Hamming distance (rjkink) in 
the A model, for p = 0.123,0.124,0.125,0.126 (curves from 
bottom to top) in case of steady state initial condition. Av- 
eraging was done over 3 x 10 5 independent samples. 

The corresponding r\ k exponent is around zero, which 
agrees with that of the PC universality class. In case 
of the survival probability we could use the conventional 
non-multi-spin coding algorithm with much less statis- 
tics. Still one can read off the same transition point with 
the value 6 ~ 0.285(8) (Fig. |), which is again in the PC 
class. 

Thus we can see the emergence of PC behaviour, which 
is in accordance with the BAWe conservation of kink- 
damage variables and the Z2 degeneration of the absorb- 
ing state arises from the fact that p c and pd coincided. 
Note, that the statistical errors are larger now than in 
the B model case, when the DS simulations were carried 
out not in the immediate neighbourhood of the critical 
point. 



The parity of the spin damage variables is not con- 
served in this case. When the simulations were started 
from random initial states we obtained the same DS tran- 
sition point as in case of the kink-damage case, but with 
neither DP nor PC universality class values. The simula- 
tions have been done both with conventional and multi- 
spin code algorithm. These resulted in the following re- 
sults for the spin damage Hamming distance as shown on 
figure Fig. 0: 
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FIG. 9. Local slopes of the Hamming distance (77) in the A 
model, for p = 0,125,0,127,0.130,0.132,0.134 (curves from 
bottom to top). The simulations were started from random 
state. Statistics over 100000 - 500000 samples. 

For the survival probability we obtained a ne arly zero 
5 exponent, as in case of kink damage; see Fig. AG: 
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FIG. 10. The same as above for the survival probability. 



These results arc quite confusing again, especially the 
exponent -q = 0.52(1), which does not belong to the 1+1 
dimensional DP or to the PC class. We can not give 
better explanation for this, that the very long transients 
prevented the healing of damages and the possibility to 
see the "true" scaling behaviour. 

Indeed, if the simulations were started from near 
steady state the results became very different. The DS 
transition point seems to coincide with the critical point 
(p c = 0.1242(5)) and we could get a spin-damage con- 
centration exponent: T) s = 0.29(2) which is close to the 
7]' = 0.285(5) of the PC scaling (Fig. |ll|). 
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FIG. 11. Local slopes of the Hamming distance (r/ s ) in the 
A model, for p = 0, 124, 0, 1245, 0.126 (curves from bottom to 
top). The simulations were started from steady state. Aver- 
aging was over 3 x 10 samples. 

The results for the survival probability (S s ) and z s co- 
incided with that of the kink-damage case, which can be 
understood by the following reasons. Although theoreti- 
cally to each spin-damage absorbing state can correspond 
two kink-damage absorbing state (by flipping all spins of 
one replica), simulations showed the kink and the spin 
damage died out always at the same time. In case of 
the spreading one can easily check that the R 2 measure- 
ments should give the same results for both kink and 



spin damage cases, because the beginning and the end 
of the perturbed region is almost the same. Indeed the 
simulations resulted in the same PC-like z (z = 1.14(1)) 
exponent in both cases (Fig 12). 
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FIG. 12. Local slopes of R 2 (t) (z) in the A model, for 
p = 0, 124,0, 1245,0.125,0.126 (curves from bottom to top). 
The simulations were started from steady state. Averaging 
was done over 5 x 10 4 independent run. 



C. Finite size scaling for both cases 

The finite size scaling simulations were performed at 
the DS transition point (pa = 0.1242) for system sizes 
L = 64, 128, ...1024. The necessary time steps to reach 
steady state were t = 10000, 20000, ... respectively. The 



results can be seen on Figure 13 
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FIG. 13. Finite size scaling results for the spin and kink 
damage for the A model. The diamonds correspond to the 
spin damage concentration, the squares to the fluctuations of 
it. The crosses correspond to the kink damage concentration, 
the "x"-s to the fluctuations of it. Triangles and stars denote 
characteristic times r of the spin and the kink damage cases. 
Averaging was done over 500 surviving samples. 

In case of the kink damage one can see regular PC-like 
scaling for the concentration: — /3fc/Vj_ = —0.5, the fluc- 
tuations of it: 7fc/^± = and for the critical dynamical 
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exponent = vn/v±_ = 1.75. In case of the spin dam- 
age we can see a constant 0.5 steady state concentration 
for all system sizes, resulting in /3 s /v± = as in case of 
the pure Glauber Ising model at T = and the NEKIM 
model at the PC transition point. In agreement with this 
and Fisher's static scaling law : 
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2f3 



(12) 



the fluctuations of it exhibit a linear scaling law ( , ~y s /b'± = 
1. Whereas the scaling of the characteristic time is de- 
scribed by the as what was found in case of the pure 
NEKIM model at the PC transition point for the spin 
variables fl27|]. 



V. THE NEKIM MODEL 

The PC universality appears in a class of non- 
equilibrium dynamic Ising models, where the kinks corre- 
sponding to '01' and TO' domain walls evolve according 
to the BAWe rules The dynamics is composed of 

the alternating application of 



a zero temperature spin flip lattice update: 



r / 7 

-(1 + 5s i ^ 1 s i+1 ) (1 - ~Si(si-i 



+ Si+l 



') 



(13) 



where 7 = tanh 2 J/kT (J denoting the coupling 
constant in the Ising Hamiltonian) , T and 5 are 
further parameters resulting in random walk, anni- 
hilation of kinks; 



and a spin-exchange lattice update: 



Wu+l 



where p ex is the probability of spin exchan 
suiting in kink — > 3 kink creation. 



(14) 



ge, re- 



The spin-flip part has been applied using two-sub-lattice 
updating, while L MC spin-exchange attempts has been 
done randomly using the outcome state of the spin- 
flip part. All these together has been counted as one 
time-step of exchange updating. (Usual MC update in 
this last step enhances the effect of p ex and leads to 
5 C = —.362(1)). In |n||27j we have started from a ran- 
dom initial state and determined the phase boundary in 
the (5, pex) plane. The phase space is composed of an 
active phase with free kinks and an absorbing, ordered 
phase without kinks provided the initial state has an even 
number of kinks. They are separated by a second order 
phase transition line of PC universality. To investigate 
the damage spreading properties we have chosen to fix 
T = 0.35, p ex = 0.3 and change S (that will play the role 
of p now). The PC critical point has been determined 
precisely f§: S c = -0.395(2). 



A. Time dependent simulations 

Time dependent simulations up to tuAX = 8192 were 
performed and we found that the DS transition point 
coincides with S c within statistical errors. The simula- 
tion results now were less sensitive whether we started 
from random initial state or from steady state. For the 
spin damage density we obtained rj s = 0.29(1) when we 
started from steady state, as in case of the A model. (If 
we started from random initial state r/ s — 0.38(2) scaling 
appeared at the 5 C point). This is shown on Fig. |lj. 




0.002 0.003 0.004 
1/t 

FIG. 14. Local slopes of the spin damage concentration 
(77), for p = 0.386, 0, 39, 0.392, 0.395, 0.4 (curves from bottom 
to top). Statistical averaging was done over 20000-40000 sam- 
ples. 

The Hamming distance measurements for kinks re- 
sulted in rj ~ in accordance with the PC universality 
class value (Fig. 15). 



0.09 
0.06 
0.03 



-0.03 



>7:: 




0.02 



FIG. 15. Local slopes of the kink damage concentration 
(jjfc), for p = 0.385, 0, 39, 0.395, 0.40.405 (curves from bottom 
to top). Statistical averaging was done over 100000 samples. 

The survival probability simulation gave the same re- 
sults for spin and kink damage cases, namely they are 
PC-like, see (Fig. |1. 
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0.015 0.02 



FIG. 16. Local slopes of the damage spin and kink survival 
probability (S), for p = 0.385,0,39,0.395,0.40.405 (curves 
from bottom to top). Statistical averaging was done over 
100000 samples. 

As in case of the A model the simulations resulted in 
the same PC-like z (z = 1.14(1)) (see Fig. [I?]) exponent 
in both cases, because the spin and kink damage regions 
start and stop approximately at the same sites. 



1.18 




0.002 



0.008 0.01 



0.004 0.006 
1/t 

FIG. 17. Local slopes of R 2 of spin damage (z), for 
p — 0.39,0.395,0.4 (curves from bottom to top). Statistical 
averaging was done over 100000 samples. 

This can be observed on Figures [l8] and [l9[ where we 
plotted the time evolution patterns of spin and kink dam- 
ages of the same run. One can see that the boundaries of 
the perturbed regions are the same, but the spin-damage 
pattern is compact, causing the higher 77^ exponent. 



B. Finite size scaling 

Finite size scaling simulations were performed at the 
DS transition point (p^ = 0.395) for system sizes L — 
64, 128, 2048. The necessary time steps to reach 
steady state were t = 7000, 14000, ...,200000 respectively. 
Figure 20. summarises the results. 




1000 1500 2000 
t 

FIG. 18. Time evolution of the spin damage in the NEKIM 
model near the DS transition point 




1000 1500 2000 
t 

FIG. 19. Time evolution of the kink damage in the NEKIM 
model near the DS transition point 



100000 
10000 
1000 
100 
10 

1 

0.1 
0.01 
0.001 
0.0001 

10 100 1000 10000 

FIG. 20. Finite size scaling results for the NEKIM model. 
The crosses correspond to kink damage concentration, the 
squares to the fluctuations of it. The stars correspond to 
spin damage concentration, diamonds to the fluctuations of 
it. The triangles and the x-s correspond to the kink and spin 
damage characteristic times (r). Averaging was done over 
500 surviving samples. The intermediate dotted line shows 
DP like scaling law for r. 



<) 



The kink concentration shows a nice scaling with 
flk/v±_ = 0.5, while the fluctuations had no pronounced 
slope on the log-log plot, suggesting -fk/v± = 0. Both 
of these values are in agreement with the PC univer- 
sality class. The spin damage concentration is constant 
(c = 0.50(1)) meaning f3 s /vx = 0.0 similarly to the A- 
model case. A linear scaling with 7 S /u± = 1 could be 
fitted for the fluctuation of it, meaning ^ s = v\_. This 
means that the exponents are in agreement with Fisher's 
scaling law (eq. [12|) again. For the critical dynamical ex- 
ponent Z = Vi\/v± PC-like scaling have been obtained for 
both the spin and kink damage cases. As the Figure |2(i| 
shows fitting can be done with L 175 (PC value) in both 
cases. The deviation from DP scaling law can clearly be 
seen on the figure too. It is reasonable to assume that the 
scaling behaviours are inherited from the pure NEKIM 
model. 



VI. CONCLUSIONS 

The damage spreading behaviour of three one- 
dimensional, non-equilibrium models, exhibiting parity 
conserving phase transition has been investigated nu- 
merically. The A SCA model of Grassberger has been 
found to be very sensitive on the initial conditions of 
the DS simulations. Acceptable results - which are in 
good agreement with that of the NEKIM model - can 
be found only when the DS simulations are started from 
steady state. 

When the DS transition point coincided (within statis- 
tical error) with the ordinary critical point of the model 
(NEKIM, A model) interesting things have happened on 
both the spin and the kink damage level. 

In case of the B model the DS transition point was 
found to be far away from the critical point and all expo- 
nents (on spin and kink damage level equally) show DP 
universality class behaviour independently of the parity 
conservation of damage variables. 

The following tables summarise the simulation results 
for the transition points and DS exponents of the mod- 
els investigated as well as DP and PC critical exponent 
estimates from refs. [1802711. 



Spin damage results 



Kink damage results 





NEKIM 


GR-A 


GR-B 


DP 


PC 


Pc 


0.395(5) 


0.1242(5) 


0.539(1) 






Pd 


0.395(5) 


0.1242(5) 


0.633(1) 






V 


0.03(3) 


0.03(3) 


0.31(2) 


0.3137(1) 


0.0000(1) 


s 


0.28(1) 


0.285(8) 


0.160(2) 


0.1596(4) 


0.285(2) 


z 


1.14(1) 


1.14(1) 




1.2660(1) 


1.141(2) 


z 


1.74(3) 


1.75(8) 


1.59(4) 


1.5798(2) 


1.750(5) 




0.500(6) 


0.48(2) 


0.27(2) 


0.2522(6) 


0.500(5) 




0.07(1) 


0.1(1) 


0.5 


0.4956(2) 


0.00(5) 


univ. 


PC 


PC 


DP 









NEKTM 


GR-A 


GR-B 


DP 


PC 


Pc 


39 c i( c >1 


1 949(^1 


VWl 1 






Pd 


395(51 


1 242(51 


633(1 1 








0.29(1) 


29(2") 


32(21 


3137(11 

U.U1U 1 111 


285(21 


s 


0.28(1) 


0.285(8) 


0.160(2) 


0.1596(4) 


0.285(2) 


z 


1.14(1) 


1.14(1) 




1.2660(1) 


1.141(2) 


z 


1.75(1) 


1.79(5) 


1.48(9) 


1.5798(2) 


1.750(5) 




0.0001(1) 


-0.001(1) 


0.26(2) 


0.2522(6) 


0.500(5) 




1.00(7) 


0.98(6) 


0.46(4) 


0.4956(2) 


0.00(5) 


univ. 


CPC 


CPC 


DP 







The rj' is the PC exponent when the BAWe process is 
started from odd number of particles. The CPC nota- 
tion denotes the compact version PC universality class 
in analogy with the CDP in case of DP universality. The 
conclusions of this paper are in agreement with previous 
works concerning PC to DP universality changes (ex- 
ample ]2~i|p7j , |3l"| ] ) . It is also in agreement with the re- 
cent DS study of Hinrichsen et. al |32|, because when 
they obtained PC class DS transition in case of the one- 
dimensional kinetic Ising model they created a mixed 
dynamics, which satisfies both the BAWe parity conser- 
vation (i.e. they followed the kinks) and the symmetric 
ground state condition (i.e. they generated passive states 
by "switching" between two dynamics that results in dou- 
ble degeneration: no-damage versus full-damage). In our 
case that has never happened that one replica would be- 
came completely reflection symmetric to the other. The 
simple dynamical rules always have driven the states to 
the same sector, where spin and kink damages died out 
simultaneously. 

One may assume this scheme to be valid for non-DS dy- 
namical transitions too, if the model has multiple absorb- 
ing states (with or without Z2 symmetry) . This hypoth- 
esis is strengthened by simulations of the pure NEKIM 
model at the PC transition point. Starting the time de- 
pendent simulations from a single seed we measured the 
usual rj and 5 exponents of the spins. In accordance 
with the DS results we obtained PC exponents again: 
77' = 0.285(5), S = 0.285(5). A detailed study of CPC 
for the sipns in the framework of NEKIM and its con- 
nection to the CDP point of Domany-Kinzel CA will be 
published soon p3| . 

These exponents satisfy the generalised hyper-scaling 
law fL§ with (3 = 0: 



£ 

3- 



2 1 + )5' + 2r}' = dz 



(15) 



where /?■ is the ultimate survival probability exponent. 
Note that the discontinuous phase transition, the com- 
pact clusters (see. Fig. |l^) and the form of the hyper- 
scaling law suggests that our case is the parity conserving 
version of the Compact Directed Percolation P9U39MQ] . 
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